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Calculations of the modal characteristics of the earth-ionosphere waveguide are presented for a 
wide range of ground conductivities. It is shown that the attenuation rate of the dominant modes, as 
a function of ground conductivity, has a maximum value which depends on frequency and various 
waveguide parameters. One of the important findings is that frequencies of the order of 10 kc/s may 
he more adversely affected by very-low conductivity than frequencies of the order of 30 kc/s. 



1. Introduction 

In most numerical treatments of VLF propagation 
in the earth-ionosphere waveguide, it is assumed 
that the ground is perfectly conducting [e.g., Al'pert, 
1960; Volland, 1961; Budden, 1962; Johler and Berry, 
1962; Galejs, 1965]. While this is surely a reasonable 
assumption for propagation over sea water, it is cer- 
tainly not valid for over-land paths [Wait, 1957]. 

It is the purpose of this paper to consider the in- 
fluence of a finitely conducting ground on the char- 
acteristics of the waveguide modes in the VLF range. 
To encompass all cases of practical interest, a fairly 
wide range of ground conductivities is considered. 
However, in order to restrict the curves to a manage- 
able number, a fixed ionospheric model is assumed 
although the reflecting height is allowed to vary some- 
what. 

2. Description of the Model 

The earth and the ionosphere are assumed to be 
spherically concentric. The earth is smooth and 
homogeneous with conductivity <j g and dielectric 
constant K (relative to free space). A convenient 
quantity to describe the characteristics of the lower 
ionosphere is the "conductivity parameter" a>, which 
is defined by [Wait, 1962; Budden, 1962] 



G0 r =<Joftlv 



(1) 



where a>o is the angular plasma frequency of the elec- 
trons and v is the effective collision frequency. A 
survey [Wait and Spies, 1964] of recent electron den- 
sity profile data shows that a reasonable daytime model 
of the lower ionosphere is specified by the relation 



Wr=(2.5X 10 5 ) exp[0(z-/?)] 



(2) 



where ft is a constant, z is the height above the earth's 
surface and h is a reference height. When z and h 



are measured in km, it would appear that j3 = 0.3 km -1 
is a representative value for the daytime ionosphere. 
At least, this is the case in the region where obliquely 
incident VLF" radio waves are reflected. At nighttime, 
it is not to be expected that this exponential form is 
adequate; however, for present purposes, this model 
is retained. 

When the effect of the earth's magnetic field must 
be considered, it is necessary to specify both the 
height profile of the plasma frequency oj () and the col- 
lision frequency v. Assuming these are exponential 
(unctions, it follows that: 



mst X exp {fi — a)z, 



: const X exp(— az) 



(3) 



(4) 



rhen 



is the inverse scale height for the collision 
frequency profile. A representative value for both 
daytime and nighttime conditions is a = 0.15 km -1 
[Wait and Spies, 1964], To display the effect of the 
earth's magnetic field in an ecomonical manner, the 
direction of propagation is chosen to be along the 
magnetic equator [Crombie, 1958]. The results are 
not quite as restrictive as they might first seem since, 
for arbitrary directions of propagation, it is known that 
the transverse component is most important for reflec- 
tion of radio waves at oblique incidence [Wait, 1962]. 
A convenient magnetic field parameter is defined by 



£1 



: ojr/^, 



(5) 



where o>r is the (angular) gyrofrequency of electrons 
in the earth's magnetic field. The quantity is nega- 
tive for propagation from the WEST TOWARDS THE 
EAST, while it is positive for propagation from the 
EAST TOWARDS THE WEST. 

For the model described above, the vertical elec- 
tric field, E, at a great circle distance, d, from a ground- 
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based vertical electric dipole source, is given by 



where E is a reference or inverse-distance field and 
W is the propagation function. Following earlier 
conventions [Wait, 1962], 



W = e~ 

(hi\) 



ITT I A 



2exp[iWC2/2]G n A M , (7) 



In the case of a flat and perfectly conducting earth, 
G n (z) may be replaced by cos (kC n z). Then, if at the 
(6) same time the upper boundary were a perfect reflector, 
it is easy to show that A n is unity. With this normali- 
zation for the excitation factor, it may well turn out 
that | A n | is sometimes greater than unity when the 
ground is finitely conducting. However, more than 
not, |A W | is somewhat less than unity when both the 
ground and the ionosphere are imperfectly reflecting. 



where X is the wavelength, A = 277-/X, G n is a height- 
gain function, C n is a complex coefficient, and A n is 
the excitation factor for mode n. The attenuation rate 
of the modes is Im(A: CJ/2) nepers per unit distance, 
while the phase velocity relative to the velocity of 
light is l + ReC*/2. 

The numerical methods used for obtaining the modal 
coefficients C„ and the complex excitation factors A /} 
have been discussed at some length in a recent tech- 
nical note [Wait and Spies, 1964], so details need not 
be given here. 

When the observer is on the surface of the earth, 
the height-gain function G n may be replaced by unity. 
When the receiving antenna is elevated one should 
refer to general expressions for the height-gain func- 
tions; however, for most cases of practical interest, 
the following expansion is satisfactory [Wait, 1962]: 



G n — 1 + i&kz — 



C\{kzf 



1 . CIA 

ka 2 



(kz)\ (8) 



where a is the radius of the earth, z is the height of 
the observer above the earth's surface, and A is the 
normalized surface impedance of the earth's surface. 
In terms of the electrical characteristics of the ground 



l€()(i> 



<j g -+- iKeoco 



1/2 



1- 



l€o(0 



erg + iKeoco 



1/2 



where co is the angular frequency and €o is the permit- 
tivity of free, space. 

When G n is near unity (i.e., small heights), it is not 
difficult to show that 



lft.Nl- 



|A|fa 
V2" 



1 (kz)\ 



where v,Jc is the phase velocity of mode n relative to 
the velocity of light. To within the same approxi- 
mation, 



phase of G n '•■ 



\L\kz 

V2" 



■a n (kz) 2 , 



where a n = k Im C|/2 is the attenuation rate in nepers 
per unit length. 

The height-gain functions are actually related to 
the excitation factor A n . In fact, the following iden- 
tity exists [Wait, 1962]: 



A n ~- 



A/2 



( h [G n (z)fdz 
Jo 



3. Discussion of Numerical Results 

Using the model of the earth-ionosphere waveguide 
as described above, extensive calculations have been 
carried out for the attenuation, phase velocity, and 
excitation factor of the modes. A fairly compre- 
hensive description of the numerical methods em- 
ployed is given in a published technical note [Wait, 
and Spies, 1964] which also contains extensive numeri- 
cal data. Here, some of the results are presented in 
graphical form with particular emphasis on the effect 
of finite ground conductivity. A number of curves 
for perfect ground conductivity are also shown to 
facilitate comparison. 

The attenuation rate (in decibels per 1000 km) as a 
function of frequency is shown in figure 1 (for mode 
number 1) for the exponential ionosphere model and 
for a perfectly conducting ground. Various reference 
or reflection heights are shown and the magnetic field 
parameter ft takes the values — 1,0, and + 1. 

It is evident from the curves in figure 1 that the 
attenuation rate for propagation from the west toward 
the east (i.e., H < 0) has the least attenuation. This 
is in conformity with much of the experimental data 
on the subject. An important feature of the curves 
in figure 1 is that the minimum attenuation for daytime 
(9) heights (i.e., h ~ 70 km) is around 16 kc/s whereas at 
night (i.e., where h ~ 90 km), the minimum attenuation 
is around 12 kc/s. These features are also quite con- 
sistent with the published experimental data [Obay 
ashi et al., 1959; Taylor, 1960; Watt and Croghan, 1964]. 

The magnitude of the excitation factor as a function 
of frequency is shown in figure 2 for the same condi- 
tions as in figure 1. The ordinate is expressed in 
decibels (dB). It is certainly evident that the excita- 
tion factor decreases with increasing frequency. In 
particular, it is seen that for the greater reflecting 
heights, the excitation of mode 1 is relatively weak. 
It is also interesting to note that the earth's magnetic 
field appears to have a very small effect on the excita- 
tion factor 

The phase velocity as a function of frequency is 
shown in figure 3 which is for the same conditions as 
figure 1. These have the expected shape of disper- 
sion curves for propagation in conventional wave- 
guides; however, one should note that the phase 
velocity v for mode 1 may be either greater or less 
than the velocity of light c. This interesting phenom- 
enon is a result of the curved nature of the waveguide. 
(12) It would appear, at least for the model employed, that 
the earth's magnetic field has an exceptionally small 
effect on the phase velocity. 
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FIGURE 1. Attenuation rate of mode I in the earth-ionosphere waveguide for perfectly conducting 

ground. 



[Here, il < corresponds to propagation from the west towards 
the west.] 
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Figure 2. The magnitude of the excitation factor of mode I in the earth-ionosphere waveguide for 

perfectly conducting ground. 
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Figure 3. Phase velocity of mode 1 in the earth-ion 

t phase of the excitation factor for mode 1 is 
i in figure 4 which, again, is for the same con- 
s as figure 1. Generally, these are positive, 
^ponding to a phase lead. Here, there appears 

some evidence of a nonreciprocal effect since 
rdinate values differ for fl = +l and ft=— 1. 
3 attenuation rate for mode 2 is shown in figure 

a perfectly conducting ground and the same 
pheric model used in figures 1-4. Again, the 
ciprocal effects are in evidence. Here, it is 
tant to point out that the attenuation rate of 

2 may be of the same order or even less than 

1 at the higher end of the VLF band. Conse- 
ly, interpretation of experimental data should 
this into account. In this paper, no further cal- 
td characteristics of mode 2 for a perfectly con- 
tig ground are shown. However, information of 
jnd may be obtained in the referenced technical 
Wait and Spies, 1964]. 

e marked influence of ground conductivity is 
ited in figure 6 where the attenuation rate for 

136 


20 24 28 8 12 16 20 24 28 
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osphere waveguide for perfectly conducting ground. 

mode 1 is plotted as a function of frequency for ground 
conductivities ranging from oo down to 0.2 millimho/m. 
In these curves, the relative dielectric constand K 
is taken to be 15. On the basis of these curves, it 
would appear that ground conductivities less than 1 
millimho/m lead to prohibitively high attenuation 
rates in the daytime (i.e., h ~ 70 km). 

The behavior of the excitation factor as a function of 
frequency, for finitely conducting ground, is illustrated 
in figure 7 for mode 1. It is rather interesting to note 
that |A| is increased by several dB as the ground con- 
ductivity changes from oo to 0.2 millimho/m. However, 
further decreases of ground conductivity, mainly at 
the higher frequencies, lead to an eventual reduction 
of the excitation factor. 

The interesting curves in figure 8 show that, for a 
height h of 90 km, the dependence on ground conduc- 
tivity is relatively weaker (e.g., compare with fig. 6). 
This is particularly the case at the higher frequencies. 
From a physical standpoint, this behavior may be at- 
tributed to the "earth-detached" character of mode 1 
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FIGURE 4. Phase of the excitation factor of mode 1 in the earth-ionosphere waveguide for perfectly conducting ground. 



in the earth-ionosphere waveguide when the reflecting 
height is of the order of 90 km or greater [Wait, 1962]. 
Thus, in this limiting situation, the energy is being 
predominantly guided along the lower edge of the iono- 
sphere without being affected appreciably by the 
presence of the earth. 

The attenuation rate and the excitation factor for 
various ground conductivities are shown in figures 9 
and 10, respectively, for h = 60 km and n=l. Simi- 
lar curves are shown in figures 11 and 12 for h = 80 km. 

A very important feature associated with finite 
ground conductivity is shown in figure 13 where the 
attenuation rate is plotted as a function of the ground 
conductivity for h = 7Q km. As indicated, the attenu- 
ation for n = 1 is seen to have a maximum value for a 
specific value of ground conductivity which itself 
depends on frequency. In general, the lower frequen- 
cies (i.e., around 10 kc/s) are characterized by increas- 
ingly high attenuation as the ground conductivity is 



decreased from 10 to 10~ 2 millimho/m. However, 
for the higher frequencies, there is a maximum below 
which the attenuation rate decreases with further 
decreases in ground conductivity. This effect is 
reminiscent of a Brewster angle phenomenon, although 
the analogy should not be carried too far since geo- 
metrical-optical concepts are of questionable validity 
here. Nevertheless, a peak in the absorption does 
suggest that the normal wave impedance in the wave- 
guide mode is nearly matched to the impedance of the 
ground. As indicated in figure 13, the effect is not 
critically dependent on the relative dielectric constant 
of the ground. 

The attenuation rate curves in figure 14 are similar 
to those in figure 13 except that now h = 90 km which 
might be more typical of nighttime conditions. The 
same maxima in the attenuation curves for n=l are 
evident. However, in this case, the ground conduc- 
tivity dependence is generally weaker, particularly 
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FIGURE 8. Attenuation rate of mode 1 for h = 90 km for various ground 
conductivities. 
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FIGURE 10. Magnitude of excitation factor of mode 1 for h = 60 km Figure 11. Attenuation rate of mode 1 for h = 80 km for various 
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at the higher frequencies. In the latter case, the 
modes have the earth-detached character which is 
consistent with this behavior. 

A number of interesting features are shown in figure 
15 where the frequency is 30 kc/s. Here, the attenu- 
ation rate and the excitation factors for modes n=l 
and n = 2 are plotted as a function of ground conduc- 
tivity for two anisotropic ionospheric models with 
h^90 km. As expected, the attenuation rate for 
n=\ shows only a weak dependence on ground con- 
ductivity. On the other hand, the ground conduc- 
tivity dependence for n = 2 is very strong and there is 
a pronounced maximum in the attenuation for a con- 
ductivity of about 0.12 millimho/m. For both modes 
1 and 2, the excitation factor shows a maximum which 
tends, to some extent, to offset the increased attenu- 
ation rates. It is also clearly evident, from figure 15, 
that the earth's magnetic field plays a relatively minor 
role here. 

The relative influences of frequency, earth's mag- 
netic field, and ground conductivity are shown rather 
clearly in figure 16 where h = 70 km. Consistently, 
the attenuation rate is least for propagation from west 
to east (i.e., fl < 0) even for relatively poor ground 
conductivities. However, it may be important to 
note that the frequency of minimum attenuation is 
significantly lowered as the ground conductivity is 
decreased from cr g =™. Of course, this effect is also 
in evidence in figure 6. 

Curves of the excitation factor, are shown in figure 
17, corresponding to the conditions of figure 16. As 
indicated, the ground conductivity has a relatively 
greater effect than that due to the earth's magnetic 
field. 



The attenuation rate for a fairly wide range of ground 
conductivities and relative dielectric constants is 
shown in figure 18. The shape of the attenuation 
versus frequency curve changes significantly for very- 
low ground conductivities. Also, in this region, the 
dependence on the relative dielectric constant of the 
ground is appreciable. 

Curves of the excitation factor, corresponding to the 
conditions of figure 18, are shown in figure 19. The 
eventual decrease of the excitation factor for very-low 
ground conductivities is certainly evident here. 

The rather interesting behavior of the phase velocity 
for low ground conductivities is seen in figure 20. In 
general, a decrease of ground conductivity is associ- 
ated with a reduced phase velocity. However, in 
certain extreme cases, the reverse seems to be the 
case. Furthermore, here, the dispersion curve has 
a negative slope which suggests that the waves are 
"backward" in the sense that phase and energy flow 
in opposing directions. However, the physical signifi- 
cance of this phenomenon is questionable since the 
associated attenuation rates are prohibitively high. 

The phase of the excitation factor, corresponding 
to the conditions in figures 18, 19, and 20, is shown in 
figure 21. It is certainly evident here that the low- 
ground conductivity leads to a substantial phase lag. 

Finally, the effect of the relative dielectric constant 
on the frequency dependence of the attenuation is 
shown in figure 22 for a relatively poor ground con- 
ductivity. The effect is really not very important, 
although it is interesting to note that the larger dielec- 
tric constants lead to lower attenuations at higher 
frequencies for n=l, while the effect is reversed at 
lower frequencies. 
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Figure 13. Attenuation rate for h=70 km as a function of ground conductivity for 

several frequencies. 



4. Concluding Remarks 

The results given in this paper should be of some use 
in making estimates of the importance of poorly con- 
ducting ground in the long-distance propagation of 
VLF radio waves. However, it should be stressed 
that the model used is highly idealized in that the 
ground is assumed to be homogeneous throughout 
the length of the path. For propagation over mixed 
paths, the situation may be much more complicated 



since discontinuities such as a land/sea boundary will 
cause conversion of modes from one order to another 
[e.g., see Wait, 1965]. However, the fact remains 
that poorly conducting portions of the actual path will 
lead to considerable weakening of the received signal 
strength. 



The authors thank William E. Garner of the U.S. 
Naval Research Laboratory, who suggested that these 
calculations be carried out. 



1367 




10 



0.05 0.1 0.2 0.5 1 2 5 

a- g GROUND CONDUCTIVITY (MILLIMHOS/METER) 
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